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ABSTRACT 

We present a measurement of the star formation properties of a uniform sample of mid-IR selected, optically 
unobscured and obscured quasars (QSOls and QS02s) in the Bootes survey region. We use a spectral energy 
distribution (SED) analysis for photometric data spanning optical to far-IR wavelengths to separate the AGN 
and host galaxy components. We find that when compared to a matched sample of QSOls, the QS02s have 
roughly twice the higher far-IR detection fractions, far-IR fluxes and infrared star formation luminosities (Tf^). 
Correspondingly, we show that the AGN obscured fraction rises from 0.3 to 0.7 between (4-40) x 10 1 Lq. We 
also find evidence associating X-ray absorption with the presence of far-IR emitting dust. Overall, these results 
are consistent with galaxy evolution models in which quasar obscuration is associated with a dust-enshrouded 
starburst galaxies. 

Subject headings: galaxies: active — quasars:general — galaxies: starburst — infrared: galaxies — X-rays: 
galaxies 


1. INTRODUCTION 

Quasars, the most luminous active galactic nuclei (AGNs 
0, have been linked to galaxies with active star formation 
(SF) ever since the discovery of their observational connec¬ 
tion to ultra-luminous infrared galaxies (ULIRGs, galaxies 
more luminous than 10 12 Lq, e.g. ISanders et al.lll988h . One 
well-studied scenario for massive galaxy evolution posits that 
gas-rich galaxy major mergers trigger both rapid supermas- 
sive black hole (SMBH) accretion and intense SF. This sce¬ 
nario associates the dust-enshrouded starburst with strong nu¬ 
clear obscuration that is later expelled by the powerful AGN, 
implying an evolutionary link between unobscured (t ype 1) 


and obscured (type 2) quasars (e.g . jDiMatteo^tal 


Hopkins et al. 2006; Gilli et al. 2007; Somerville et al. 


2005 


20081 


Treister et al.1120091k 

On the other hand, the “unification model” of AGN as¬ 
cribes obscuration of AGNs to different lines of sight through 
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a dusty “torus” surround ing the SMBH (e.g. lUrrv & Padovanil 
1 1 995llAritonuccil 1 1 9931) . This model predicts no difference 
in host galaxy properties between obscured and unobscured 
AGNs. To date, it is still a matter of debate whether the 
obscuration in luminous quasars can be explained solely by 
the orientation-based unification model or if it is also en¬ 
hanced due to dust on larger scales throughout the host galaxy. 
Several studies have shown results supporting a scenario de¬ 
parting from the unification model, such as the enhanced 
SF activity (e.g. ICanalizo & Stocktonll200lllPaue et al.ll2004b 

iHiner et al1l2009t iBrusa et al.ll2008|: IShan & Chen 2012]) and 

the more disturbed structure (e.g. iLacv et al.ll2007bl) of the 

host galaxies of dust-obscured quasars when compared to un¬ 
obscured quasars. Clustering of different types of AGNs have 
also shown that obscured A GNs are more strongly clustered 
than unobscured AGNs (e.g.lHickox et al.l2011tlDonoso et al.1 
1201 4t iDiPomneo et al.l 120141) . However, other studies have 


found no significant difference between obscured and unob¬ 
scured AG N populations in the i r morphological and SF prop- 


erties (e.g. JSturm et al. 2006; Zakamska et al. 

20061 2008; 

Mainieri et al. 2011; Schawinski et al. 2012b 

Merloni et all 

2013), and the host galaxy star formation rate (SFR) does 
not distinguish X-ray selected AGNs with different obscur- 

ing column densities (e.g 

Rovilos et al.1120121 Rosario et al.l 

2012: Merloni et al. 20l| 

. 


Nonetheless, the measurements of SF properties for quasar 
host galaxies still suffer from selection biases that are often 
different among various quasar populations. In particular, op¬ 
tical and X-ray selected quasar samples might have different 
completeness in obscured and unobscured sources due to the 
attenuation of optical and X-ray radiation by dust and gas. In 
fact, some studies have suggested that optical surveys might 
miss ~ 50% of the AGN population due to both the obscura¬ 
tion a nd host galaxy contamination in low-luminosity AGN s 
(e.g. IGoulding & Alexander! I2009t iGoulding et all 120101) . 
While current X-ray observations probing photons with en¬ 
ergy at ~ 10 keV are believed to be less affected by moderate 
levels of obscuring materials (Ah < 10 24 cm -2 ), a significant 
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fraction of X-ray AGNs (~ 20-50%, e.g.lDon lev et al.l|2005 


Guainazzi et al.l 20051:iPark et al.ll2010t [A lexander et al.ll201 1 

Georganto poulos et al.1120 1 3t IWilkes et al.1 120 1 3l) might still 


be missed due to obscuration that is Compton-thick (Ah 
10 24 cm -2 ). Although very hard X-ray photons can penetrate 
Compton-thick obscuration, the previous high energy X-ray 
surveys are limited to local sources due to the shallow flux 
limits (e.g. Swift/BA T, iBurlon et al.l 1201 ll and INTEGRAL, 
Sazon ov et al.l 120121) . While the NuSTAR (lHarrison et al.1 


20131) mission opened a new window of high energy X-ray 


up to 80 keV, the recent studies that used NuSTAR to observe 
heavily obscured AGNs have found that obscuration is still 
a non-n egligible effect even for X-ray photon s at such high 
energy (IStern et al.ll2014t ILansburv et alJl2014l) . 

In contrast, mid-IR observations of the reprocessed emis¬ 
sion from the obscuring dust can detect heavily obscured 
AGNs. A number of studies have shown that large pop¬ 
ulations of AGNs can be selected using the power-law 

sn gths (e.g 

. 120071: Donley et al 


2004 

; Stern et al. 

2012 

; Assef et al. 


2005t IHickox et al 


2010b). Although these mid-IR color 


selection criteria cannot avoid some star forming galaxy 


interlopers and mig ht miss AGNs ac c reting at lower ac- 


cretion rates (e.g. Hickox et al.l 2009 

: | Donley et al. 120121 

Mateos et al.l2013t Hainline et al. 2014 

Chung et al. 2014!) or 


AGNs with complicated silicate features in their mid-IR SED 
([Kirkpatrick et al.ll2013l) . they are effective in selecting both 
obscured and unobscured AGNs with similar completeness in 
mid-IR wavelengths. Recent studies have shown that mid-IR 
selected AGNs can also be separated into obscured and unob¬ 
scured populations u sing a simple optical to mid-IR color se- 
lection criterion (e.g. IHickox et al .1120071 iDonoso et al.ll2014t 

DiPo mpeo et alj|2014l) . which can easily be explained by the 

different level of extinction in the optical emission from the 
nucleus (e.g. IHickox et al.ll2007l H07 hereafter). Therefore, 
a direct comparison of the host galaxy SF properties of ob¬ 
scured and unobscured sources in a mid-IR selected quasar 
sample can provide insights to the origin of obscurations in 
rapidly accreting SMBHs. 

In this work, we adopt the mid-IR selected quasar sam¬ 
ple from IHickox et ali ( 201 ll Hll hereafter), which is com¬ 
prised of unobscured and obscured quasars with similar distri¬ 
butions in quasar properties (e.g. redshift and quasar luminos¬ 
ity), thus making it an excellent sample to study the connec¬ 
tion between host galaxy properties and quasar obscuration. 
The quasar sample studied in this work consists of 546 unob¬ 
scured quasars (QSOls) and 345 obscured quasars (QS02s) 
selected using Spitzer mid-IR observations in the Bootes sur¬ 
vey region. We utilize the optical spectroscopy from the AGN 
and Galaxy Evolution Survey (AGES. lKochanek et al.11201 2l) 
and th e XBootes Chandra X-ray observations (Murray et al l 
2005), along with the 250 ^m data f rom the Spectral and 
Photometric Imaging Receiver (SPIRE, Grif fin et al.l20 10) on 
board the Herschel Space Observatory. With the inclusion 
of the far-IR photometry in which AGN contributions have 


been shown to be much smaller than that of the host galaxy 

(e.g. 

Netzer et al. 20071; Facy et 

al. 2007a; Kirkpatrick et al. 

2012 

; Mullanev et al. 2012llChen et al.ll2013t Del Moro et al. 

2013 

; Drouartet al. 2013), we can use the wealth of multi- 


wavelength observations in the Bootes field to obtain robust 
measurements of SFR for luminous quasars and test if they 
are associated with obscuration. 


We also take advantage of this sample to study the cor¬ 


relation between star formation and AGN accretion for the 
mid-IR QSOs. In galaxy evolution models which suggest 
that the concurrent growth of AGN and host galaxy occurs 
during the dust enshrouded galaxy merger phase, different 
types of quasars represent different stages of galaxy evolu¬ 
tion. To date, observational studies on the connection be¬ 
tween SFR and AGN accretion rate have not yet shown defini¬ 
tive conclusions. Positive, strong correlations have been 
found in studies of optically selected type I quasars (e.g. 
Serjeant & Hatziminaogloul 120091) and type II quasars (e.g. 
Netzerl2009l) : while different conclusions have been found in 
X-ray selected quasars at higher redshift (e.g.lSilverman et al.l 
l2009t iMullanev et al.ll2012t iRosario et alJl2012l) . The differ- 
ence between observed correlations could be driven by var¬ 
ious limita tions and obse r vational constraints such as sam¬ 
ple size (Page et al.l 120121 l Harrison et al.l 120 121) or AGN in - 
trinsic variability (e.g. IChen et al.ll2013fc IHickox et al.ll2014h . 
In addition, it is also important to note that AGNs selected 
with different crit eria are hosted by very different host galaxy 
populations (e.g. IHickox et al.ll2009t iGriffith & SternfeoIOt 
iGoulding et al.l 20141) in which the galaxy and SMBH might 
follow different evolutionary paths. The sample in this work 
consists of a large number of luminous QSOs including both 
the unobscured and heavily obscured populations, thus bi¬ 
ases due to small number statistics and the exclusion of dust- 
enshrouded AGNs are small. The bias in the observed SFR- 
AGN accretion rate relation in AGN host galaxies due to 
short-term AGN accretion rate stochasticity mav_ also be less 
significant in luminous quasars IHickox et al.ll2014l) . There¬ 
fore, the SFR-AGN accretion rate correlation for the mid-IR 
selected quasar populations may shed light on the origin of 
the concurrent growth of SMBH and galaxy. 

This paper is organized as follows: in §2 we describe the 
multi-wavelength data and the properties of the quasar sam¬ 
ple. In §4, we discuss the SED fitting procedures that we used 
to disentangle the AGN and the host galaxy contributions. In 
§3, we explore details of the observed far-IR properties of 
QSOls and QS02s. A comparison of the SF luminosity be¬ 
tween QSOls and QS02s is laid out at §5; The x-ray proper¬ 
ties for the QS02s are discussed in §6 and the AGN obscured 
fraction as a function of SF luminosity is discussed in §7. A 
discussion and a summary are given in §8. Throughout the pa¬ 
per, we use the Vega magnitude system and assume a ACDM 
cosmology with f l,„ = 0.3, = 0.7 and Ho = 70 km s _I . 

2. THE QUASAR SAMPLE 

In this section, we discuss the data from the Bootes survey 
region as well as the quasar selection and classification criteria 
we adopted. 


2.1. Data 

The sample studied in this work comes from the 9 deg 2 
Bootes survey region covered by the NOAO Deep Wide-Field 
Survey (NDWFS, iJannuzi & Devi U999I) . Bootes is unique 
among extragalactic surveys because of its large area and the 
excellent multi-wavelength coverage from space- and ground- 
based telescopes, which make possible statistical study of the 
rare luminous AGNs. 

In this work, we use the multiwavelength photometry cat¬ 
alog from Brown et al. (private communication), which the 
same the one used in lChung et al.l (120141) . This catalog cov¬ 
ers the optical to mid-IR bands including the NDWFS opti¬ 
cal observations in Bw, R. I bands, near-IR photometry from 
the NOAO NEWFIRM survey (J, H and Ks, iGonzalez et al.l 
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Fig. 1.— (a) IRAC color-color diagram showing the selection of the quasar samples using the criteria of IStem et alJ <20051) . The gray-scale shows the density 
of sources detected at > 5cr significance in all four bands in SDWFS data. Blue stars and red circles show the QSOl and QS02 samples, respectively. The 
I Stem et all <20051) color-color selection region is shown by the dashed line. In addition, the QSOs with far-IR detections (FIR-QSO) are enclosed with orange 
circles. Some of the mid-IR QSOs fall out of the selection wedge due to the updated IRAC photometry and aperture corrections (see §2.3). (b) Illustration of 
the optical-IR color-selection criteria for dividing the IR-selected QSO sample into unobscured (QSOl) and obscured (QSQ2) subsample s. Shown is observed 
/? —[4.5] color versus bolometric luminosity, calculated as described in §4. Contours show the distribution for all the lHickox et al.l <20071) IR-selected quasars, 
while blue stars and red circles show the QSOl and QS02 subsamples at 0.7 < z < 1.8 used in HI 1 and this analysis as described in §2. The distribution in the 
R — [4.5] color and Lagn ar e also shown as histograms in the side panels. QSOls are shown as the blue solid lines, QS02s are shown as the red dashed lines and 
the FIR-QSOs are shown as the orange dash-dotted lines. The contours and color histograms show that a simple cut in optical-IR color clearly separates the QSO 
samples into two populations. 


2010b and the Spitzer Deep Wide Field Survey (SDWFS, 
Ashby et al.l2009l) of the 4 bands of mid-IR observations from 
the Spitzer Infrared Array Camera (IRAC) at 3.6,4.5,5.8 and 
8 pm. In addition, the mid-IR photometry from the Spitzer 
Multi-band Imaging Photometer (MIPS) at 24 pm is also in¬ 
cluded (IRS GTO team, J. Houck (PI), and M. Rieke). In this 
catalog, the optical to near-IR photometry was extracted from 
a matched aperture for each band. For Bw, R, I, H and Ks 
bands, the photometry was measured from images smoothed 
to a common point spread function (PSF) with a 1735 full 
width half maximum (FWHM); while for / band, the photom¬ 
etry was measured from images smoothed to a common PSF 
with a 1 ."60 FWHM. To ensure consistency in the photome¬ 
try, we use 6"aperture photometry from optical bands through 
the IRAC bands in the mid-IR to account for the large IRAC 
beam size. For the MIPS data, we use PSF photometry in¬ 
stead of aperture photometry due to the still larger beam size 
of MIPS. The 5<t flux limits of the optical to near-IR broad¬ 
band photometry are 25.2, 23.9, 22.9, 21.1, 20.1, 18.9 (Vega 
magnitudes) for the Bw, R, I J, H, Ks bands, respectively. For 
the mid-IR wavelengths, the 5cr flux limits are 6.4, 8.8, 51, 50 
and 170 pJy for the 3.6,4.5,5.8,8.0 and 24 pm bands, respec¬ 
tively. A extensive description of the multiband photometry 
extraction can be found in lBrown et al . (1200 7). 

We also make use of the far-IR observations from the 
publicly available Herscliel Multi-tiered Extragalactic Survey 
(HerMES. lOliver et al.ll201 2l) . We re-reduced and mosaiced 
the Bootes SPIRE observations (lAlberts et alJ 120131) . which 
include a deep ~ 2 deg 2 inner region near the center of the 
field and a shallower ~ 8.5 deg 2 outer region. We specifically 
focused on removing striping, astrometry offsets, and glitches 


missed by the standard pipeline reduction. We also convolved 
the raw maps with a matched filter (see iChapin et allf20111) . 
which aided in source extraction by lowering the overall noise 
and de-blending sources. From this, we generated a matched 
filter catalog with a 5er detection threshold. In this catalog, we 
consider SPIRE sources with fluxes larger than 20 mJy as un¬ 
ambiguously detected. Completeness simulations show that 
these catalogs are 95% complete in the inner region and 69% 
complete in the outer regions above a flux of 20 mJy. We 
also find minimal flux boosting for low SNR sources above 
this flux cutoff (see S2.2 in lAlberts et al.ll20l3l for the details 
of the completeness simulation). We match the positions of 
the SPIRE catalog to the /-band positions with a matching 
radius of 5". We tested the rate of spurious matches by offset¬ 
ting the SPIRE source positions by F in a random direction 
and matching the randomly shifted catalog to the Bootes cata¬ 
log. We found that with a radius of 5", our matching between 
the SPIRE and Bootes catalog only yielded < 2% spurious 
matches. 

In addition, we also match the positions of the lBrown et al.l 
(120071) catalog to the publicly available Wide-fwld Infrared 
Survey Explorer (WISE) All-sky catalog ( iWright et al.ll2010l) 
and obtain the profile-fit photometry magnitudes in the Wl, 
W2, W3 and W4 bands (3.4,4.6,12 and 22 pm). 

As a complementary measurement of the AGN accretion 
rate and absorption by gas, we utilize X-ray data from the 
XBootes survey, whi ch is a mosaic of 126 short (5 ks) Chan¬ 
dra ACIS-I images dMurrav et al.l 12005) iKenter et~akl 120051) 
covering the entire NDWFS. XBootes contains 3,293 X-ray 
point sources with at least four counts in the AGES survey 
region. The conversion factors from count rates (in counts 
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s 1 ) to flux (in erg s 1 ) for the XBootes are 6.0 x 10 12 erg s 1 
count -1 in the 0.5-2 keV band and 1.9 x 10 11 erg s -1 count -1 
in the 2-7 keV band, which are derived for a 5ks on-axis ob¬ 
servation and assuming a_canonical unabsorbed AGN X-ray 
spectrum (see §3.3 of lKenter et al.ll2005h . 

For this study, we use spectroscopic redshifts ( spec-zs ) from 
AGES when possible. For the sources without a spec-z , we 
adopt the photometric redshifts ( photo-zs ) calculated using 
techniques combining artificial neural network and template 
fitting algorithms (iBrodwin et al.l [20061) . The uncertainty of 
this set of photo-z is a = 0.06(1 +z) for galaxies and a = 
0.12(1 +z) for AGNs. We note that the uncertainty for AGNs 
are dominated by the difference between photo-zs and spec- 
zs for the type I AGNs. The uncertainties in the photo-zs for 
QS02s are difficult to estimate accurately as only 9% of the 
QS02s have spec-z measurements. For these 32 QS02s, the 
uncertainty of the photo-zs is o = 0.06(1 +z), which is consis¬ 
tent with the uncertainty for the galaxy population. Recently 
spectroscopic follow-up of 35 W7.S7i-selected, optically ob¬ 
scured quasars by Hainlin e et al.l (120141) al so found that photo- 
zs derived from template-based algorithm (lAssef et al.l 20 10a) 
have a similar accuracy. This is not surprising since the opti¬ 
cal SED for the heavily obscured AGNs is dominated by the 
host galaxy which has more spectral features and should allow 
for a more accurate photo-z measurement. 

An upper limit for the QS02 photo-z uncertainties can be 
estimated based on different approaches. H07 obtained an un¬ 
certainty of a z = 0.25(1 + z) by comparing the IBrodwin et al.l 
( 2006 ) photo-zs to the photo-zs estimated by fitting 3 different 
galaxy templates to the Bw, R, I photometry. Since the accu¬ 
racy of the photo-zs based on fitting 3 galaxy template to 3 op¬ 
tical photometric observations is limited, the o z = 0.25(1 + z.) 
uncertainty is a very conservative upper limit. For this study, 
we adopt the photo-z uncertainty upper limit of <j z = 0.25( 1 +z.) 
from H07. However, this is a very conservative estimate 
since the photo-z uncertainties for the QS02s with host galaxy 
dominated optical SEDs are likely to be much smaller. A full 
discussion of the photo-z uncertainties is given in §8. 

2.2. AGN identification and classification 

The AGNs in this work is drawn from the quasar sample of 
HI 1, which is a subset of quasars from the H07 mid-IR AGN 
sample. The H07 AGNs were identified based on the IRAC 
color-color selection criterion of IStern et al.l d2005l) . For the 
redshift range of the H07 sample (0.7 < g < 3.0), the IRAC fil¬ 
ters probe wavelengths at which the characteristic power-law 
continuum from the reprocessed AGN radiation starts to dom¬ 
inate. At these wavelengths, the light from old stellar popula¬ 
tions is characterized by a Rayleigh-Jeans tail of a blackbody 
radiation with temperature higher than 2500K, which peaks 
at wavelengths different than both the AGN accretion disk 
and the reprocessed AGN emission. In addition, dust emis¬ 
sion at near-sublimation temperature heated by the AGN is 
significantly stronger than that heated by massive stars, thus 
AGNs can be easily identified with their mid-IR col or (e.g. 
lLacv et akll2004tIStern et al.l2005tlDonlev et al.ll2012l) . These 
reprocessed photons at mid-IR wavelengths are less affected 
by obscuration than the optical and the X-ray photons due to 
their much smaller absorption cross section. Therefore, even 
heavily obscured AGNs can be identified using mid-IR selec¬ 
tion criteria. 

An important feature of the H07 catalog is that the 
AGNs were selected from the IRAC Shallow Survey (ISS, 
lEisenhardt et al.l[20041) . which has shallower flux limits than 


the later SDWFS survey. While the IStern et ahl (120051) AGN 
selection can be contaminated by SF galaxies in deep mid-IR 
surveys (e.g. lDonlev et al. 2012), the con tamination is negli¬ 
gible at the shallow flux limits of ISS (see lHickox et aklliZOfWl 
lAssef et al.ll2010ail20TlT) . 

H07 also showed that the distribution of the optical ( R ) 
band to mid-IR (4.5pm) colors is bimodal for the luminous 
mid-IR quasars. This bimodality can be easily explained 
by the difference between the SEDs of unobscured and ob¬ 
scured quasars. For mid-IR selected quasars, the AGN dom¬ 
inates the mid-IR wavelengths, while for obscured quasars, 
the nuclear emission at optical wavelengths is heavily ab¬ 
sorbed. Thus the SED in the R band for obscured quasars is 
similar to those of normal galaxies (e.g. iPolletta et al.l 120061: 
iHickox et al.l 120071 and Hickox et al. 2014, in preparation). 
H07 has shown that the mid-IR selected quasars can be sepa¬ 
rated into two distinct populations of quasars with an empiri¬ 
cal color cut at R- [4.5] = 6.1. Most (~ 80%) of the quasars 
in H07 with R- [4.5] <6.1 are spectroscopically confirmed as 
type 1 quasars, therefore in this work we refer to these unob¬ 
scured quasars withR-[4.5] < 6.1 as QSOls. As for quasars 
with R— [4.5] > 6.1, due to the limited depth of AGES, only 
~ 5% have spectroscopic measurements and can be classified 
as type 2 quasars spectroscopically. However, H07 have ex¬ 
tensively studied the quasars with R— [4.5] >6.1 and have 
shown that these quasars are consistent with bright, obscured 
X-ray AGNs with Ah > 10 22 cm -2 and with AGN bolometric 
luminosities Lagn > 10 45 erg s -1 . Therefore, we refer to the 
obscured quasars with R— [4.5] >6.1 as QS02s. 

This empirical classification based on the quasar mid-IR to 
optical color has been shown to be broadly consistent with 
both the spectroscopically classified sources dDonoso et al.l 
120141) and X-ray hardness ratio (defined as HR = (H-S)/(H+ 
S ) where H and S are the photon counts in the 2-7 keV 
ban d and 0.5-2 keV band, respectively) cl assification crite¬ 
ria (IHickox et all 120071: lUsman et al.l l2014h . Since we lack 
spectroscopic confirmation for the majority of the QS02s, it 
is possible that the QS02s are different than the optical type 
2 quasars. Recently. lLacv et ahl (12013l) have also studied the 
optical and near-IR _spectroscopy of mid-IR AGNs selected 
using the lDonlev et al.l (120121) IRAC color criteria, and found 
that a significant fraction (~ 33%) of the AGNs with evident 
mid-IR power-law continuum have no significant emission 
lines consistent with the optical emission line diagnostics (e.g. 
iBaldwin et al.l 1 19811) . suggesting that these mid-IR obscured 
AGNs are more deeply obscured than the optically selected 
type 2 AGNs^ This has also been suggested by the recent 
study of lHainline et al.l (120141) . 

2.3. The final sample 

To minimize the uncertainty in the measurements of quasar 
and host galaxy properties due to the scattering in the photo- 
zs, we focus on the Hll quasar sample, which is a subsam¬ 
ple from H07. The Hll quasar sample focuses on the red- 
shift range 0.7 < z < 1.8 and limits the QSOls to those with 
broad optical emission lines and robust spec-z measurements 
from AGES. For the QSOs in this redshift range, all of the 
sources have I < 18, and the spectroscopic QSOl sample is 
highly complete. For the QS02s in this redshift range, only 32 
(~ 9%) of the QS02s have spec-z mea surements. For the rest 
of QS02s, we adopt the photo-zs from lBrodwin et ahl (120061) . 
This subset of QSOs have been used in the study of clustering 
properties of QSOls and QS02s in Hll, which showed that 
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in the redshift range 0.7 < z < 1.8, the sample is highly com¬ 
plete at the shallow flux limits of ISS. The mid-IR selection of 
this particular sample is shown to suffer less than < 20% con¬ 
tamination from star-forming galaxies dHickox et all 12011 ). 
Since the purpose of this work is to compare the star forma¬ 
tion properties of QSOls and QS02s, it is important to make 
sure that any difference in SF properties is not driven by the 
presence of starburst galaxy interlopers. Also, even for bright 
mid-IR quasars the starburst contribution in the MIR may be 
non-negligible. Therefore, we rely on SED decompositions to 
disentangle the AGN and starburst components. The details 
of the SED decomposition are laid out at §3. Our SED fitting 
results show that all of the Hll QSOs have a non-negligible 
(> 20%) AGN component. Based on the SED fitting results, 
~ 3% of the HI 1 QSOs would have Lagn less than the QSO 
criterion (Lagn > 10 45 erg s -1 ). We therefore limit our focus 
to the sample of 546 QSOls and 345 QS02s which satisfies 
the QSO luminosity criterion and comprise ~ 97% of the Hll 
sample. 

As a demonstration of the sample used in this work, we 
show the distribution of the IRAC [3.6] - [4.5] to [5.8] - [8.0] 
colors and the R — [4.5] to Lagn distributions of the quasar 
sample in Fig. Q] We note that in Fig. [TJi, a small number of 
the H07 sources have [3.6]-4.5 and [5.8]-[8.0] colors that 
lie outside t he IStern etal.l (120051) AGN selection wedge. The 
IStern et all (120051) AGN selection wedge and the H07 sam¬ 
ple were defined using the original ISS catalog in which the 
mid-IR photometry was derived using the standard aperture 
correction, which is different than the PSF profile fitting cor¬ 
rections for individual sources in the [Brown et al. ( 2007 ) cat¬ 
alog. Moreover, the updated IRAC photometry comes from 
SDWFS, which is 2 times deeper than the ISS. Thus, it is not 
surprising that a small fraction (~4%) of the original H07 
sample does not meet the origin al ISter n et ~l. (1200 5) criterion. 
However, the QSO sample in Fig.[TJi shows a tight [5.8]-[8.0] 
to [3.6]- [4.5] distribution, which is not seen in the similar 
Fig. la of Hll which use the original ISS photometry. This 
suggests that the majority of our sample does have power- 
law like mid-IR SED with the high-precision photometry of 
SDWFS (e.g. see Fig. 22 in H07). For the purpose of com¬ 
pleteness we do not exclude the sources outside of the updated 
AGN selection wedge, since these sources might still have a 
non-negligible AGN component which can be identified by 
the SED fits. However, we note the exclusion of this small 
fraction of sources has little effect on the average properties 
of the QSOs. We also show the redshift and flux distributions 
of the sample in Fig. [2] which shows that the redshift distri¬ 
butions of the QSOls and QS02s are similar. 

3. SED DECOMPOSITION 

Although it is now well-established that far-IR fluxes 
at wavelengths longer than 60 pm in AGN host galaxies 
are mostly dominated by the emission related to SF (e.g . 

iNetzer et al.ll2007b iMullanev et al.ll20llt iRosario et al.ll2012D . 

caution is still required when estimating the SF properties of 
powerful quasars. Here we use SED decompositions to ensure 
that we have reliable estimates of the intrinsic AGN luminos¬ 
ity (Lagn) and star formation luminosity (L!jF~T.) 

To properly disentangle the emission from stars and AGN 
accretion, we fit all of the available broad band photometry in 
our datasets with SED templates. We combine several empir- 

13 Defined as the integrated 8 — 1000 fim luminosity of the host galaxy 
component only. 



Fig. 2.— Top: The 250 gm distributions for QSOls and QS02s with direct 
SPIRE detections. QSOls are shown as the blue stars. The 5 far-IR detected 
QS02s with spec-z measurements are shown as the filled circles. The rest 
of the far-IR detected QS02s with photo-zs are shown as open circles. Bot¬ 
tom'. The redshift distributions of QSOls and QS02s are shown in blue and 
red solid lines. Also shown are the redshift distributions of the QSOls and 
QS02s detected at 250 p m. The histograms for far-IR detected QSOls and 
QS02s are shown as the blue and red dash-dotted lines and filled in orange. 
This plot shows that the redshift distributions of QSOls and QS02s are sim¬ 
ilar and QS02s have more 250 //m detections than QSOls. 

ically derived templates to create AGN and host galaxy tem¬ 
plates covering the wavelength range from optical to far-IR in 
our data (see §3.1). 


3.1. SED templates 

Our SED fitting procedures utilize the four empirical SED 
templates spanning 0.03 - 30 pro described in lAssef et ail 
(1201 OaL A10 hereafter). A10 have shown that the non¬ 
negative combinations of three different galaxy templates and 
a single AGN template, with the addition of extinction to 
the AGN component only, can robustly describe the opti¬ 
cal to mid-IR SEDs of a wide variety of AGNs selected us¬ 
ing Spitzer IRAC color (lA ssef et al .l l2010al) or WISE color 
(lAssef et al.ll2010bt IChung et al .1120 1 4f) . However, these tem¬ 
plates do not include the far-IR wavelengths that are crucial 
for this work. 

To extend the A10 AGN template to the far-IR, we create 
ad hoc AGN templates by replacing the hot-dust component 
of the A10 AGN template with the average infrared quasar 
template from lNetzer et al.l d2007|) and the three infrared AGN 
templatesjfor low-, mean-, and high-luminosity AGNs) from 
IMullanev et al.l (1201 ll) . which cover a wide range of AGNs 
with different mid-IR to far-IR properties. This is based on 
the assumption that the SED of the hot accretion disk around 
the SMBH is the same for all the radiatively efficient AGNs, 
and that the differences are due to variable extinction and dif¬ 
ferent spectral shapes in the mid- and far-IR wavelengths due 
to different distributions of dust. To account for the extinc¬ 
tion at the AGN templates, we use the iDraine d2003l) extinc¬ 
tion law with Ry =3.1, which mainly attenuates the SED at 
A < 30 pm and also produces the silicate absorption features 
at mid-IR wavelengths that are common among AGN. The ex¬ 
tinction strength is treated as a free parameter with a range of 
0 < A v <48. 

For the host galaxy templates, we consider two different 
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components: the contribution from the stellar population of 
the galaxy, which accounts for the optical to near-IR emis¬ 
sion; and a starburst component, which represents the mid- 
to far-IR dust emission from re-processed stellar light. For 
the stellar population component, we follow the approach 
described in A10 by assuming that the stellar population in 
any galaxy is comprised of the non-negative combination of 
the three empirical galaxy templates with populations of star- 
burst (Im), continuou s star-forming (She) and old stars (el¬ 
liptical), respectively ( lAssef et al.ll2008l . : 20l0al) . Unlike the 
elliptical template, the Sbc and Im templates both contain hot 
dust components in addition to the stellar population SEDs. 
Since the dust component does not extend to the far-IR wave¬ 
lengths we require and will be taken into account in the star- 
burst templates which we will choose later on, we replace the 
dust emission (A > 4.9 /im) in the Sbc and Im templates with 
the SED identical to the elliptical galaxy to create empirical 
stellar population templates with no dust emission. At these 
wavelengths the stellar SEDs are dominated by low mass stars 
similar to the elliptical template. For our mid-IR selected 
quasar sample, the star light contribution is negligible at the 
wavelengths where hot dust emission dominates. 

For the starburst component, we use the 105 starburst tem¬ 
plates from lCharv & Elba?] (120011) an d the 64 starbur s t tem¬ 
plat es fromlDale & Heloul (120021) . The lCharv & Elba?! (120011) 
and iDale & Heloul (120021) starburst templates cover a wide 
range of SEDs for various prototypical local star forming 
galaxies, with Lir in the range between 10 s - 10 135 L 0 . For 
SF galaxies at z > 1, these templates have been shown to be 
reliable when used to estimate SF-related Lir with monochro¬ 
matic SPIRE ob servations (e.g. |E baz et al .11201 It) . However, 
as pointed out by Kirkpatrick et al. (2012j), the starburst galax¬ 
ies at higher redshift have different dust temperatures and 
SED shapes. Although the effect of the different dust temper¬ 
ature on the estimation of L^ is small when the wavelengths 
of the SPIRE bands probe the peak of the cold dust emission, 
it is still important to include high redshift starburst templates 
to accommodate the possibly different SED shapes of high 
redshift starbursts. Therefore, we also include th e z ~ 1 and 
z ~ 2 a verage starburst SED templates from Kirk patrick et al J 
(120121) in our SED fitting analysis. We thus adopt a total of 
171 starburst templates in our SED fitting analysis. 

Given these SED templates described above, we fit the ob¬ 
served photometry with an iterative \' 2 minimization algo¬ 
rithm (Levenberg-Marquardt) to minimize the function. 

^filters 

x 2 =£ 

,=° (i) 

^obs,i — 6^star,i — & CXt(ii (B — V))^AGN,i — £*^starburst,i 

^ ' 

Here F 0 b Sj i is the observed flux in the ;th band and a, is its 
uncertainty. F ga | axy j, Lagnj and -Fstarburstd are the fluxes of the 
stellar (optical), AGN (optical to far-IR) and starburst (far-IR) 
templates for filter i. The function ext gives the extinction 
of the AGN flux in band i given the color excess E(B-V). We 
optimize a.b.c and E(B-V) to minimize the y 2 . 

3.2. Results 

While uniform NUV to mid-IR photometric coverage is 
available for all of the sources in our sample, only ~ 18% 
of them have a SPIRE 250 /jm detection with a flux larger 


than the 20 mJy. Since our sample is selected to have AGN- 
dominated mid-IR SED, the far-IR photometry is essential to 
constrain the average SF properties. For the sources with di¬ 
rect far-IR detections, we fit their multiwavelength photome¬ 
try with the SED templates described in the previous section. 
For the 731 sources without direct far-IR detections, we use a 
stacking analysis to constrain their average far-IR flux. 

We stack the QSOs without direct far-IR detections in bins 
of R— [4.5] (see §4.2) and bins of Lagn (derived in HI 1 from 
interpolations between the IRAC photometry and the MIPS 
24 /rm photometry). The uncertainty of the stacked flux is 
determined by bootstrap resampling. We created 10,000 ran¬ 
dom samples by drawing objects from the original samples 
with replacement until the number of objects in each random 
sample is the same as the number in the original sample. The 
uncertainty of the stacked flux is the variation in the stacked 
fluxes of the random samples. The details of the far-IR slack¬ 
ing analysis can be found in §3.1.1 of Alberts et al.1 ( 2013 )). 
The results of the stacking analysis are given in Table 1. 

To estimate the average SF luminosity of the far-IR non- 
detected QSOs, we assume that their 250 /rm fluxes and 
250 /rm flux uncertainties are equal to the averages found 
from the stacking analysis for the sources with similar values 
of both R— [4.5] and Lagn- We find that for the majority of the 
SPIRE-non detected sample, the best-fitting SED reproduces 
the observed average far-IR flux well. 

From SED fitting results with far-IR fluxes stacked in bins 
if R - [4.5] and Lagn, we find that ~ 85% of the quasars 
have a prominent AGN hot dust component (AGN component 
> 50%) at mid-IR wavelengths covered by the IRAC bands. 
All of the quasars have a non-negligible (> 20%) AGN com¬ 
ponent, confirming that the Hll sample consists of power¬ 
ful quasars. However, the AGN dominates the quasar SED 
even at 24 /im for more than a half of our sample. So the SF 
properties of these powerful quasars is primarily constrained 
by their far-IR photometry. For individual sources with only 
the stacked SPIRE flux to anchor the starburst component, 
this approach of using stacked far-IR flux can lead to inac¬ 
curate Lj^. However, the main purpose of this work is not 
to determine the L^ for individual QSOs. Instead, we take 
this approach to constrain the AGN contamination in the av¬ 
erage far-IR SF luminosity with the well-fitted mid-IR AGN 
SED, which can robustly determine the average Lj^ for dif¬ 
ferent QSO populations. We note that when we compare the 
stacked fluxes from either the R— [4.5] bin or the Lagn bin 
with the best-fitting AGN component, we find that the aver¬ 
age AGN contribution at 250 /rm is only 5%; and none of the 
QSOs in our sample have an AGN component that contributes 
more than 40% at 250 /rm. Since the average intrinsic SEDs 
for AGNs are expected to fall more rapidly in more luminous 
systems (iMullanev et al.llTol ll) . a starburst component is es¬ 
sential to reproduce to average flux at 250 /rm. 

With the best-fitting SEDs, we calculate the total infrared 
luminosity for each QSO (L[r ) by integrating the best-fitting 
SEDs from 8- 1000 /rm. For each source in our sample, we 
also calculate the Lir of the starburst component by integrat¬ 
ing the host galaxy component of the best-fitting SED over 
the same wavelength range (Lj^ hereafter) [3- We find that 
for our sample, the average ratio between L 3 ,' and L \r is 

14 For clarification, we use Lir as a general terminology for any integrated 
8-1000 /rm luminosity. Lj^* is the integrated Lir of the best-fitting SED and 
Lis the Lir of the host galaxy component only. 
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Fig. 3.— Examples of the best-fitting SEDs (solid line) in the rest-frame of each source. These sources are fitted using an AGN component with a IDrainel 
(200 3!) extinction law (dotted curve), an empirical stellar component (dash-dotted) and empirical starburst templates (dashed). We find that almost all sources 
have an SED dominated by the starburst component in the far-IR, see §3 for details. 


~ 55%, indicative of a non-negligible AGN contribution to 

in t hese mid-IR luminous qu asars. This is also consistent 
with the [Kirkpatrick et all d2012b results of the study combin¬ 
ing Spitzer IRS spectra and multi-band Herschel photometry 
for z ~ 1 and z ~ 2 ULIRGs. However, the far-IR flux probed 
by the SPIRE 250 pm band is dominated (> 70%) by the star- 
burst component for ~ 90% of our sample. This result con¬ 
firms that quasar SEDs are still dominated by a starburst com- 
ponent at rest-frame wavelengths longer than 100 pm (e.g . 
iNetzer et alJl2007t iMullanev et alJl201lt iRosario et alJl2012t 

iDel Moro et alJ20T3h . Therefore, the average L^ of powerful 
mid-IR quasars can still be robustly measured with the inclu¬ 
sion of far-IR photometry. 

We next derive the bolometric AGN luminosity (Lagn) by 
directly integrating the de-absorbed AGN component. We 
compare the Lagn from our SED fitting with the Lagn derived 
in HI 1, and found that our Lagn is higher than the HI 1 L A gn 
by an average of ~ 0.1 dex. Since the derivation of Lagn in 
HI 1 did not correct for the weak but non-negligible dust ob¬ 
scuration and SF galaxy contamination in the mid-IR, it is not 
surprising that our Lagn is different. However, the difference 
is small due to the fact that most of the quasars in this sam¬ 
ple are dominated by the AGN and the dust obscuration only 
weakly attenuates the SED at mid-IR wavelengths. We show 
the redshift distributions of L\'^ and Lagn in Fig. [4] 

4. FAR-IR OBSERVATIONS AND QUASAR OBSCURATION 

Several studies have pointed out that even for AGNs with 
quasar-like luminosites, their far-IR SEDs are often dom- 
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Fig. 4.— The redshift distributions of the (8—1000 fim) Lj^ 31 and the AGN 
bolometric luminosity derived from SED fitting and far-IR stacking for non- 
detected sources described in §3.2. The symbols and lines represent the same 
subsets of objects shown in Fig. 1, where blue stars and red circles are QSOls 
and QS02s, respectively. Individually far-IR detected sources (FIR-QSO) are 
enclosed with orange circles. 


inate d by the cold dust emission heated by young stars 
(e.g. Netzer et all 120071 lLacv et all [2007 at iKirkpatrick et al.1 
12012 ; IMullanev et al.ll2012h . Although alternative cases have 
been reported i n some recent studies (e.g. iHiner et aD 120091: 
iDai et al.1l2012l) . the cold temperature of the far-IR emission 
still requires the AGN-heated dust to reside at a large distance 
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R-[4.5] (Vega) 

Fig. 5.— The 250 fin l detection fraction (/ 250 ) as a function of the R — 
[4.5] color. The /250 for all QSOls is shown as the blue star, and the /250 
for all QS02s (weighted to have a redshift and AGN bolometric luminosity 
distribution similar to that of QSOls) is shown as the red circles. This figure 
shows that QS02s are 2.4 times more likely to have a far-IR detection. 

from the central SMBH ( Sanders et al.l 1 988 ). In addition, 
powerful molecular outflows might also produce strong far- 
IR emission (e.g. iSun et all 120141) . However, the number of 
AGNs and starburst galaxies with confirmed warm molecular 
outflows from high-resolution far-IR or sub-mm observations 
is limited. Nonetheless, strong far-IR emission for quasars 
implies the existence of dust extending to a large distance 
well beyond the putative obscuring torus, and a comparison 
of far-IR properties between QSOls and QS02s could de¬ 
termine that whether the quasar obscuration is related to the 
large-scale dust. 

4.1. SPIRE detection fraction 

We begin with a very simple test by measuring the far-IR 
detection fraction above a flux limit of 20 mJy in the Her- 
schel SPIRE 250 //m filter (/250 hereafter) for QSOls and 
QS02s separately. In our survey region, the shallowness of 
the SPIRE observation implies that any quasars in our sam¬ 
ple with a far-IR detection are at least as bright as luminous 
infrared galaxies (LIRGs, which are defined as galaxies with 
Lir(8 — 1000 gm) > 1 O U L 0 ) at the same redshifts. At first 
glance, we find that among the 546 QSOls, only 68 (12%) of 
them are detected in the SPIRE 250 gm filter, while 102 of 
the 345 (29%) QS02s have far-IR detections. 

We next use R— [4.5] as a rough proxy of the obscuration 
strength on the nuclear emission and study whether /250 is re¬ 
lated to the extinction at the R band. The primary goal of this 
analysis is to study whether the star formation properties of 
the quasar-hosting galaxies are related to the observed obscu¬ 
ration. Therefore, it is important to make sure the QSOls and 
QS02s are matched in key properties, i.e. redshift and AGN 
bolometric luminosity (Lagn, estimated from the SED fitting 
analysis discussed in §3). We divide our sample into four bins 
of R— [4.5] and calculate /250 for each bin by applying statis¬ 
tical weights to each quasar so that the samples are matched 
in redshift and Lagn in each bin. We find that the weighted 
/250 increases rapidly with R - [4.5] for QSOls. For QS02s, 
/250 is only weakly correlated with R — [4.5], but the /250 of 
QS02s is higher than that of QSOls by a factor of 2.7 ± 0.2. 
We show the /250 to R- [4.5] relation in Fig. [5] 

4.2. The average SPIRE 250 gun fluxes 


Due to the flux limit of the SPIRE catalog, /250 only reflects 
the fraction of QSOs hosted by galaxies with very luminous 
far-IR emission. To extend this analysis, we estimate the av¬ 
erage 250 /itn flux in bins of R— [4.5] to test if the far-IR flux 
of our quasar sample also evolves with R— [4.5]. The measure 
the average 250 gm flux (S250) f° r all QSOs, we first esti¬ 
mate the average 250 /itn flux for the sources without direct 
SPIRE observations using the stacking analysis discussed in 
§3.2. Next, the stacked fluxes is combine with the sources 
individually detected in far-IR to calculate the average S 250 
for all QSOs. We show the average 250 gm flux (S250) in the 
same R— [4.5] bins used for Figure [5] with results given in 
Table 1. 

From the stacking analysis, we find that for QSOs without 
direct far-IR detections, there is no significant dependence of 
the average S 250 on R - [4.5]. However, the mean .So 50 of the 
individually detected far-IR sources shows a R-[4.5] depen¬ 
dence similar to that of / 250 . Therefore, driven by the increas¬ 
ing fraction of far-IR luminous sources, the average .S 250 of 
the entire QSO sample is correlated with R - [4.5]. 

This result shows that for powerful mid-IR quasars with 
more dust attenuation at optical wavelengths, the average far- 
IR emission is stronger. This implies that at least part of the 
obscuration seen in the obscured QSOs might be due to the 
far-IR emitting dust. 

5. THE AVERAGE Lf^ OF QSO IS AND QS02S 

In §3, we found that QS02s have higher /250 and S 250 than 
QSOls. In this section, we examine whether the observed dif¬ 
ference in the average far-IR fluxes of QSOls and QS02s can 
be associated with the difference in the star formation proper¬ 
ties of their host galaxies, as predicted by the galaxy evolution 
models with a connection between quasar obscuration and star 
formation. 

We first calculate the median star formation luminosity 
by averaging the derived from SED fitting (§4). We 
find that for QSOls, the mean is lo 45 - 29±003 erg s 1 
(1O U - 71 L 0 ); and for QS02s, the mean is 10 45 - 59± ° 04 
erg s -1 (10 12 01 Lq). We find that similar to /250 and the av¬ 
erage S 250 , the median for QS02s is significantly higher 
than that of QSOls by 0.30 dex. The higher Lf £ in QS02s 
confirms the results from §3 that QS02s are associated with 
host galaxies with more star-forming cold dust which might 
be obscuring the nuclear emission. 

5.1. Lfl vs. R —[4.5] 

To further study the relation between AGN obscuration and 
, we divide our QSO sample into bins of R - [4.5] as a 
proxy for AGN obscuration. To be consistent with §3.1 and 
avoid the uncertainties by the lack of spec-zs of QS02s, we 
opt to use R— [4.5] instead of the E(B-V) derived from SED 
fitting as the proxy for nuclear obscuration, since the uncer¬ 
tainties of the photo-z limit the accuracy of our ability to mea¬ 
sure E(B-V). We calculate the median AGN contribution at 
observed frame R and [4.5] bands using the best-fitting SEDs 
in each bin. We find that towards redder R - [4.5] colors, 
the AGN contributes 89%,79%,29%, 12% at the R band and 
74%,59%,43%,60% at the [4.5] band. This suggests that the 
empirical R— [4.5] color is indeed a good proxy of AGN ob¬ 
scuration at optical wavelengths, and that the QSO classifica¬ 
tion based on the R— [4.5] color is reliable. We then calculate 
the average in bins of R— [4.5]. We again estimate 

the uncertainty in ) using bootstrap resampling. 
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TABLE 1 

Results of Far-IR Stacking Analysis 


Average /250 in bins of R-4.5 

*-[4.5] (Vega) 

4.7 

5.4 

6.6 

7.5 

S 25 o(mJy) (ND) 

6.58 ±0.73 

7.93 ±0.59 

7.39 ± 0.75 

8.57 ±1.00 

S 2 50(mJy) (All) 

8.04 ±0.84 

11.4 ±0.77 

16.4 ± 1.50 

18.0 ±0.21 

(z) (all) 

1.31 

1.20 

1.19 

1.31 

Average /250 in bins of Lagn 

log (Lagn) [ergs '] 

45.2 

45.6 

46.1 

46.5 

S 25 o(QS01XmJy) 

3.76 ±0.82 

5.49 ±0.70 

5.21 ±0.75 

7.70 ±0.95 

S 25 o(QS02)(mJy) 

3.56 ±0.94 

8.84 ±1.40 

5.18±0.66 

8.44 ±2.70 

z(average) 

0.98 

1.17 

1.36 

1.49 


Note. — Results from the far-IR stacking analysis described in §4.2. The top half of the table shows the results for sources binned in R — [4.5], The first 
row shows the results for the far-IR non-detected (ND) sources only, and the second row shows the result for the entire QSO sample. The average redshift for all 
QSOs in each bin is also given. The bottom half of the table shows the results for the entire QSOl and QS02 samples binned in their AGN bolometric luminosity. 
The average redshift for each R— [4.5] and Lagn bin is also listed; we note that there is only a < 0.02 average redshift difference between the QSOl s and QS02s. 

46.0 

45.8 
|j> 45.6 

a j 
c/3 ^ 

m 45.4 
o 

45.2 

45.0 

4 5 6 7 8 45.5 46.0 46.5 47.0 

R-[4.5] (Vega) log L AGN [erg/s] 

Fig. 6.— (a) Left : The relationship between and R-[4.5] for the mid-IR selected QSOs. Far-IR detected QSOs are shown by the orange, open symbols and 
the stacking results for the far-IR non-detected QSOs are shown by the gray, open symbols. In both cases, there is only marginal increase in L§j with R— [4.5]. 
However, driven by the much higher far-IR detection fraction in QS02s, the average for QS02s is higher than the average Lj^ for QSOl. (b) Right : the 
average for QSOs in bins of Lagn- The results for QSO Is are shown as the filled, blue stars; while the results for QS02s are shown as the filled red circles. 
For the entire QSO population, we also show their average Lj^ in four bins of Lagn as the open, green squares. For comparison, we show the lHickox et alj (20141) 
model evaluated at 0.7 <z< 1.8 as the shaded region. We note that the average redshift in each bin increases with Lagn! and the relationship between (Lj^) and 
Lagn is consistent with the redshift evolution for main-sequence SF galaxies with M* = 10 u Mq. We show the redshift evolution of MS galaxies as the dashed 
line. The L^ for all QSOs coincides with the evolution of MS galaxies, suggesting a connection between QSO host galaxies and MS star-forming galaxies (see 
§5.2). 
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We show the results in Fig. [ 6 ji. Even though the dif¬ 
ference between QSOls and QS02s is large, (Lf^) does not 
show the strong dependence on R - [4.5] seen for /250 and the 
average S 250 (Fig- Eh ) within QSOls and QS02s. Also, the 
/-i'll, difference between the far-IR detected QSOl s and QS02s 
and the far-IR non-detected QSOls and QS02s are only 0.14 
dex and 0.17 dex, which are not as significant as the 0.30 dex 
difference between all QSOls and QS02s. Thus the higher 
(Lj^) of the QS02 sample is mainly driven by its higher frac¬ 
tion of far-IR luminous SF galaxies. 

5.2. /-['ll, vs. Lagn 

To explore the connections between AGN and host galaxy 
growth rates in different QSO populations, we separately di¬ 
vide the QSOls and QS02s into bins of Lagn, and mea¬ 
sure the (Lj’r) for each bin. The uncertainties are again 
estimated by bootstrapping as discussed in §5.1. We plot 
the Lagn-L^ relations for QSOls and QS02s in Fig. Eh- 
We find that both the QSOl sample and the QS02 sample 
have a positive Lj^-Lagn correlation with similarly shallow 


slopes, with log oc 0.25 logLAGN for QSOls and logLj^ oc 
0.271ogL AG N for QS02s. However, the (Lj^) for QS02s is 
higher than that of QSOls in each Lagn bin by 0.28,0.25,0.32 
and 0.28 dex. 

For comparison, we also show the Lagn — L®^ model from 
iHickox et al.l(l2014l) which assumes a direct, linear correlation 
between the average Lagn and L*^ while the observed Lagn is 
modulated by rapid variability. We calculate the lHickox et all 
d2014l) relation spanning the 0.7 < z < 1.8 redshift range in 
Fig. [6] as the shaded region. We find that for all QS02s be¬ 
sides those in the most luminous Lagn bin, our result is con¬ 
sistent with the H14 model, and the QSOls have Lj^ slightly 
lower than the model predictions throughout the entire Lagn 
range. However, for both types of QSOs, the slope of the Lj^- 
Lagn relation appears to be shallower than the prediction of 
the simple model if we take the different average redshift in 
each Lagn bins into account. 

As was shown by the results in Table 1, the average red- 
shifts for more luminous QSOs are higher. Although our far- 
IR stacking approach can reliably measure the average far-IR 
fluxes and hence the average Lj’ | | , it is important to note that 
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the difference in the average redshift between the lowest and 
the highest Lagn bins is 0.5, and it is possible that the ob¬ 
served correlation between and Lagn is partially driven 
by the cosmic evolution of star formation and AGN accretion. 
Even though recent studies of the evolution of cosmic infrared 
luminosity density have shown that for the redshift range of 
the four Lagn bins in Fig. [6J3, the cosmic infrared luminosity 
is only weakly correlated with redshift (pnt ex (1+z) -0 3±tu for 
1.1 < z < 2.85. lGruppioni et~aPl2013l) . there is still a substan¬ 
tial evolution in the specific SFR (SFR per unit stellar mass ) 
of “main sequence (MS)” SF galaxies (e.g. iElbaz et al .112(1 1 it) . 

A simple estimate of the effect of MS galaxy SFR red¬ 
shift evolution can be made by calculating the average Lj^ 
for MS galaxies with M* = 10 u Mq at the redshift range of 
our QSO sample. From the theoretical framework of dark 
matter halo abundance matching methods (e.g. lBehroozi et akl 
120131) . 10 u Mq is the typical galaxy stellar mass hosted by 
dark matter haloes of mass M/ m / c = 10 131 [/r’M Q ], which 
is the halo mass r eported in recent c lustering studies of 
mid-IR QSOs (e.g. iHickox et all 1201 lb iDonoso et akl l2014t 
iDiPomneo et al ll20l 4) . At the redshift of each QSO in our 
sample, the Lj^' for a 10 h Mq MS galaxy can therefore be 
evaluated using the redshift-dependent SFR-M* relation from 
(W hitaker et al.l (120121) and a Kennicutt relation, SFR = 1.09 x 
L^?/L ? (lKennicutl T998. modified for a Chabrier initial mass 
function). We find that the average Lj^ for MS galaxies eval¬ 
uated using this approach is consistent with the average Ljr 
in each Lagn bin. This implies that for mid-IR quasars in 
this redshift range, the Ljr and instantaneous Lagn might not 
be directly connected but still follow similar redshift evolu¬ 
tion. This suggests that a common physical parameter (e.g. a 
common gas supply) might be driving the evolution of both 
SMBHs and galaxies. This is also consistent with the recent 
study of iRosario et al J (12013!!) . which showed that the Lf^- 
Lagn correlation for broad emission line QSOs is consistent 
with a scenario in which quasars are hosted by normal star¬ 
forming galaxies. 

For mid-IR selected AGNs, contamination from star¬ 
forming galaxy interlopers is also an issue that might intro¬ 
duce biases to the measurements of AGN and star formation 
luminosities. While our SED fitting results show that 95% 
of the sources in our sample are indeed powerful QSOs with 
AGN-dominated mid-IR SED, it is still very important to ver¬ 
ify that the observed Lj^-Lagn relation in Fig.[6]is not caused 
by the incompleteness of our AGN selection criterion. In 
particular, a very powerful SF galaxy can harbor a heavily 
obscured AGN and still have starburst-like mid-IR SED. In 
suc h case, the I RAC color of the SF galaxy would fall out of 
thelStern et akl (12 0051) selection wedge (e.g. [Kirkpatrick et al.1 
l2012tlAssef et al. ll2010bH Chung et all 120141) . and the sample 
selected with the lStern et al.1 (120051) wedge would show a bi¬ 
ased Lj^-Lagn relation. 

By happenstance, the mid-IR color of typical star-forming 
gal axy templates occupies the lower-left corner of the 
IStern et a k| (|2005| ) wedge at the redshift range of this sam¬ 
ple (e.g. iDonley et al.l 12012L Fig. 17). Therefore, any AGN 
contribution in addition to the SF galaxy templates would eas¬ 
ily promote a SF galaxy into the AGN identification wedge. 
This is particularly true for our sample of luminous quasars. 
We test this effect by first normalizing the 6 /im luminosity 
for AGN templates we used in the SED fitting described in 
§4 to L A gn = 10 45 erg s -1 to meet the minimum Lagn crite¬ 


rion of our QSO sample. We next combine these normalized 
AGN templates with two archetypical starburst galaxy tem¬ 
plates, M82 and Arp 220 from lPolletta et al . (200 7). We find 
that to have a star forming galaxy falling out of the AGN 
identification wedge while having an Lagn ~ 10 45 erg s -1 , 
its Lj^ must exceed 10 47 erg s -1 . This suggests that an SF 
galaxy must be a hyper-luminous infrared galaxies (HyLIRG, 
Lir > 1 0 i3 Lq) to be able to hide an underlying quasar-like 
AGN component. According to the far-IR luminosity func¬ 
tion of lGruppioni et al.l d201 31) , there would only be less than 
2 HyLIRGs in the redshift range of our sample given the size 
of the Bootes survey region. Because of the large number of 
sources in the lowest L A gn bin of our sample, the inclusion of 
the 2 additional QSOs hidden in HyLIRGs would only affect 
the average by 0.15 dex, which is still within the 2 a range 
of the original Lsf- 

We can also estimate biases in Lf^ due to the sources not 
selected in the AGN identification wedge by examining the 
SEDs of the sources not classified as mid-IR QSOs. In the 
redshift range of our sam ple, there are 77 SPIRE-detected 
sources outside the lStern et al.) ( 2005 ) wedge. When we fit the 
SEDs of these sources we find that most of these powerful SF 
galaxies have no significant AGN contribution in the mid-IR. 
For the sources without a mid-IR AGN component, we as¬ 
sume that AGN contributes < 10% to the mid-IR monochro¬ 
matic luminosity at 6 /im as an upper limit. Only 4 of the 77 
have Lagn satisfying the Lagn > 10 45 erg s -1 criterion. If we 
add these four sources into the lowest Lagn bin, the average 
would only increase by 0.05 dex. Therefore, we conclude 
that the Lsf-Lagn correlation observed in Fig. [6] is not bi¬ 
ased by the exclusion of heavily obscured quasars hidden in 
powerful SF galaxies. 

6. X-RAY PROPERTIES OF MID-IR QSOS 

In §5, we have shown that for mid-IR selected QSOs, the 
AGN obscuration at optical wavelengths and the AGN mid-IR 
luminosity can both be connected to the star formation of their 
host galaxies. Here we use an alternative AGN accretion rate 
indicator, the X-ray emission, to study the interplay between 
the far-IR emitting dust and the X-ray emission. 

We first count the number of detections in the far-IR and the 
X-ray observations. We note that the counts of X-ray detec¬ 
tions can be affected by the varying sensitivity of the Chan¬ 
dra observations across the field (e.g. Mendez et al.112013 ). 
However, the XBootes observations are relatively uniform in 
depth (~ 4-8 x 10 -15 erg s -1 . iKenter et al.1 2 0(')5l i and this 
effect should be equivalent for both QSO types. The results 
of the detection fractions are summarized in Table 2. First, 
QSOls have an X-ray detection fraction of 65%, which is 
much higher than the 12% far-IR detection fraction for QSO 1 s 
(see §3). For QS02s, the X-ray detection fraction is only 
34% while the far-IR detection fraction is 27%. We also find 
that the presence of far-IR emission is associated with lower 
X-ray detection fractions in both QSOls and QS02s. For 
QSOls, 91% of the X-ray detected AGNs are not detected 
in the far-IR; while for QS02s, 73% of the X-ray AGNs have 
no detectable far-IR emission. These number counts imply 
that large scale dust might also play an important role in the 
absorption of X-rays. 

Since a large fraction (~ 76%) of the far-IR detected QS02s 
have no direct X-ray detection, these far-IR bright QS02s 
might not be included in an X-ray selected AGN sample. 
However, since both the far-IR and the X-ray observations 
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TABLE 2 

Number of X-ray detected sources 



Nx 

Nnx 

fx 

(z) 

log {Hum ) [ergs -1 ] 

(log U*)(L Q ] 

(GrHergs- 1 ] 

(Av) 

QSOl (ALL) 

356 

190 

0.65 

1.24 

46.20 

11.74 

44.17 ±0.03 

0.32 

QSOl (FIR) 

30 

35 

0.46 

1.26 

46.38 

12.21 

43.93 ±0.09 

0.50 

QSOl (No FIR) 

326 

155 

0.67 

1.24 

46.17 

11.62 

44.19 ±0.03 

0.28 

QS02 (ALL) 

119 

235 

0.34 

1.26 

46.25 

12.06 

43.78 ±0.05 

2.62 

QS02 (FIR) 

23 

74 

0.24 

1.24 

46.27 

12.34 

43.66 ±0.10 

2.85 

QS02(No FIR) 

87 

161 

0.35 

1.26 

46.24 

11.80 

43.82 ±0.05 

2.59 


Note. — The number counts and properties of different subsamples of QSOs classified based on the X-ray and far-IR detections. The properties listed here are 
the number of sources with X-ray detections (Nx), the number of source without X-ray detections (Nnx), the X-ray detection fraction (fx )■ the average redshift 
( (z)). the average AGN bolometric luminosity ((log Lagn)» see §4). the rest-frame 2-10 keV X-ray luminosity and Ay (see §4 and §5.3). 


in our sample are relatively shallow, the detection fractions 
only reflect the incidence of bright starburst and bright X- 
ray AGNs. To compensate for the shallow flux limit of the 
XBootes survey, we use an X-ray stacking analysis to estimate 
the average X-ray luminosity for the sources without direct 
X-ray detections. We defined the stacked X-ray counts as the 
average number of background-subtracted photons detected 
within the 90% point-spread function (PSF) energy encircled 
radius at 1.5 keV, rgo, where rg 0 = 1"+ 1O"(0/1O') 2 - Here 6 is 
the off-axis angle from the Chandra optical axi^O We adopt 
3.0 and 5.0 counts s _l deg -2 for the diffuse X-ray background 
at 0.5-2 keV and 2-7 keV bands, which closely matches the 
average background in annuli after excluding detected X-ray 
sources. We refer to §5.1.1 ofH07 and lChen et all (120131) for 
details about X-ray stacking analysis. 

We performed the analysis on the subsamples listed in Ta¬ 
ble 2. To calculate the average X-ray luminosity for each sub¬ 
sample of QSOs, we first calculate the rest-frame 2-10 keV 
Lx for each X-ray detected quasar from their observed 0.5- 
7 keV luminosity with a ^-correction using the spectral in¬ 
dex derived from the flux ratio between the soft (0.5-2 keV) 
and the hard (2-7 keV) bands. In detail, the X-ray spectrum 
was assumed to be a simple power-law and the hardness ratio 
was used to determine the power-law index (T) using PIMMS 
and the ACIS Cycle 4 on-axis response function. For the X- 
ray detected QSOls, the average index is T = 1.83, which is 
typical for unabsorbed AGN X-ray spectra. For the X-ray de¬ 
tected QS02s, the average index of 1= 1.24 is consistent with 
AGNs with moderate absorbed spectra (Hickox et al.]f2007 '). 
With the r derived from the observed hardness ratio, we k- 
corrected the observed-frame 0.5-7 keV luminosity to the 
rest-frame 2-10 keV, without de-absorbing the X-ray spectra. 
For the X-ray non-detected sources, we calculate their aver¬ 
age rest-frame 2-10 keV Lx by /^-correcting the stacked 0.5-7 
keV luminosity derived from the stacked X-ray flux using the 
median redshift and the hardness ratio of stacked X-ray pho¬ 
tons. The general X-ray properties of QSOls and QS02s have 
been discussed at length in §5 of [Hickox et al. ll ([200 7;). Com¬ 
bining the Lx of the individually X-ray detected sources with 
the average Lx for the non-detected sources from stacking, we 
obtain the average X-ray luminosities in Table 2. These X- 
ray luminosities are the observed values and are not corrected 
for absorption. We also list the median extinction Ay of the 
AGN template at optical wavelengths from the SED fitting re¬ 
sults (see §4). For heavily obscured AGNs, the SED at optical 
wavelengths can be dominated by the host galaxy, and the av¬ 
erage Ay may be a lower limit. However, we do find a much 
larger value of A y for far-IR detected QSOls than that for far- 

15 Chandra Proposers’s Observatory Guide (POG), available at 
http://cxc.harvard.edu/proposer/POG 


IR non-detected QSOls, which implies that the presence of 
star-forming dust attenuates both the optical and X-ray AGN 
emission. 

We can also use the total numbers of stacked hard-band 
(2-7 keV) and soft-band (0.5-2 keV) photons to obtain the 
hardness ratio for different QSO subsamples. We find that for 
QSOls, the HR remains a constant HR ~ -0.42 for both the 
far-IR detected and far-IR non-detected subsamples. While 
for QS02s, the HR for both the far-IR detected and far-IR 
non-detected subsamples are also very similar (HR ~ -0.12). 
The large difference in the HR between QSOls and QS02s 
is consistent with the typical HR values for type I and type 
II X-ray AGNs, which has been demonstrated in H07. How¬ 
ever, we here point out that the average hardness ratio is in¬ 
sensitive to the presence of far-IR emitting dust. The uni¬ 
form HRs for far-IR detected and far-IR non-detected QS02s 
might simply be due to the wide redshift distribution from 
which the stacked sources are drawn from, since the uncer¬ 
tainties of the measured HRs corresponds to more than 0.3 
dex in N H (2-5 x 10 22 cm -2 , e.g. Fig. 15 in H07) for the red¬ 
shift range of our sample. However, the hardness ratios for the 
far-IR detected and far-IR non-detected QSOls are both con¬ 
sistent with little or no neutral gas absorption, which cannot 
explain the large Lx difference between far-IR QSO 1 s and far- 
IR non-detected QSOls. Therefore, the constant HR regard¬ 
less of the presence of strong far-IR emissions might imply 
that the lower average Lx and the higher Ay in far-IR detected 
subsamples are due to gas of large column density that ob¬ 
scures X-ray photons in both hard-band and soft-band. We 
note that large columns of ionized gas residing within lOpc 
of the vicinity of the X-ray emitting corona might also ab¬ 
sorb both the hard-band and soft-band photons without al¬ 
tering the observed HR. As d emonstrated in various s tudies 
of X-ray spectral models (e.g. Ceballos & Barcons 1 9%h. the 
observed hardness ratio is a much weaker function of the col¬ 
umn density of the ionized warm absorber when compared to 
a neutral, cold absorber. This specific type of absorption has 
been observed in QSOs hosted by extremely luminous sub¬ 
millimeter galaxies (e.g. [Page et al.ll201lT) . similar to the lu¬ 
minous SF galaxies hosting QS02s in this work. However, 
whether the same type of warm absorber exists in our far-IR 
detected QSOS cannot be confirmed without high-resolution 
X-ray spectroscopy. 

We summarize the results of our stacking analysis in Fig. [7J 
where we show the ratio between L^/m and (Lx) for dif¬ 
ferent subsamples of QSOs. Since Lx and Lmir are both 
excellent tracers of SMBH accretion and Lmir is relatively 
insensitive to obscuration, the Lx to Lmir ratio has been 
used to track down the elusive population of Compton-thick 
AGNs and to study the intrinsic AGN accretion rate (e.g . 
iLutz et all l2004t lAlexander et all 120081: iGandhi et all 120091 
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Fig. 7.— A summary of the Lx/Le^m for different QSO populations. For 
comparison, we show the Lx /As/mi relationship for type I X-ray AGNs using 
Eq. 2 from IFiore et alJ (20091) evaluated at Z^m = lO^ ^erg s -1 as the dashed 
line. We also sho w the Lx/L^um relatio n for local AGNs with Compton- 
thick obscurations dAlexander et alJl2008D as the dotted line. We show that 
the presence of far-IR emitting dust attenuates the observed Lx in comparison 
to Lg^rn an d increases the AGN extinction in the optical wavelengths. We 
note that far-IR luminous QS02s might not be detected in wide-held X-ray 
surveys due to the heavy obscuration that can be associated with the dust- 
enshrouded host galaxies. 

Eckart et al.H2010t IGoulding et alj|201 It iRovilos et all 120131: 
Lansburv etal.ll2()14l IStern et al. 1120 141). 

In Fig. H we show that for both QSO Is and QS02s, the 
existence of a far-IR detection is connected with smaller 
L x /ratios by at least 0.2 dex when compared to the 
far-IR non-detected subsamples. This indirectly supports a 
link between far-IR emitting dust and the absorption of X- 
ray photons. For comparison, we also show the luminosity- 
dependent L(,j ,m —Lx relation bv lFiore et al.l (l2009t) evaluated 
at the average of our sample = 10 449 erg s” 1 . We 
find that the L x /L^pm for QSOls without direct far-IR de¬ 
tections is similar to the iFiore et all ( 2009 ) relation. We note 
that for local active galaxies, an almost linear G,/;m - Lx re¬ 
lation has been reported by several studies (e. g. iLutz et aD 
l2004t iMaiolino et al.l 120071 iGandhi et al.l [20091) , but the in¬ 
ferred Lx/Ls^m value at L^pra = 10 44,9 erg s -1 would be -0.37 
if we adopt Eq. 2 from iGandhi et all (120091) . which is much 
higher than the results for our far-IR non-detected QSOls and 
the IFiore et akl (2009) relation. The X-ray hardness ratio of 
the QSOl sample is consistent with little or no X-ray absorp¬ 
tion, suggesting that the different Lx/L^pra ratios for QSOs 
and the X-ray selected AGNs (IFiore et alj|2009t) observed are 
not simply due to the difference in gas absorption. The differ¬ 
ence in the intrinsic mid-IR to X-ray spectral shape between 
local AGNs and quasars is beyond the scope of this work but 
will be addressed in a follow-up work by Chen et al. (2015, 
in preparation). 

For AGNs with Compton-thick absorption, the observed Lx 
can be attenuated by a factor of more than 15 when com¬ 
pared toJheinfrinsicLg^m-Lx relation (e.g. I Alexander et akl 
l2008t IGoulding et al.l2010l) . We find that the Gy/Ga/m of our 
far-IR detected QS02s is consistent with the Alexand er et al.l 
( 2008 ) relation derived from attenuating the X-ray luminosi¬ 
ties for the L x /G//m relation of local AGNs (ILutz et al.l 
12004] ) with Compton-thick material. However, recent obser¬ 
vations have also discovered several luminous AGNs with ex¬ 
tremely weak X-ray emissi ons even in the ultra-hard energy 
bands probed by NuSTAR dLuo et al.ll2013l:lTeng et al1l2014t 
Lansb urv et al.ll2014t IStern et al .1120141) ."Therefore, it is still 


unclear whether the low L x /Gi/im ratio for the QS02s is due 
to Compton-thick obscuration or intrinsic X-ray weakness. 
Nonetheless, the low L x /G,//m for far-IR detected QS02s is 
another line of evidence indicating that X-ray AGN selection 
methods can miss heavily obscured mid-IR QSOs hosted by 
galaxies with active SF activity. We caution that all AGN se¬ 
lection criteria suffer from various different selection biases, 
especially the incompleteness due to obscuration and the con¬ 
tamination from SF galaxy interlopers. We reiterate that SF 
galaxy contamination with our mid-IR color selection crite¬ 
ria is relatively mild for the shallow mid-IR flux limits of 
our sample. Our result demonstrates the advantage of mid-IR 
AGN selection criteria in detecting deeply obscured quasars 
hosted by powerful SF galaxies. 

7. OBSCURED FRACTION AND Lf[ 

In §5 and §6 we have shown that mid-IR selected QS02s 
reside in galaxies with more active SF and that their AGNs are 
more obscured at both optical and X-ray wavelengths. While 
these results support the connection between AGN obscura¬ 
tion and star formation in the QSO host galaxies, we next 
directly test the connection between quasar obscuration and 
galaxy SF with the measurement of the “obscured fraction” 
as a function of L^. We show this result in Fig. [8] For the 
far-IR detected QSOs, where can be more accurately esti¬ 
mated, we separate the sources into two bins of with equal 
number of sources in each bin. For the QSOs with esti¬ 
mated from stacked far-IR fluxes, the L^ suffer from larger 
uncertainties. Therefore, we do not group these QSOs based 
on their Lf^ but only evaluate obscured fractions for the entire 
populations. We find that the obscured fraction of the far-IR 
detected QSOs is much higher than that of the rest of QSOs, 
and the obscured fraction increases monotonically with re¬ 
spect to Lf^ as / 0 bs oc 0.39 x logLj^. This agrees with a sce¬ 
nario which connects the obscuration in powerful quasars to 
the host galaxy star formation, as the AGN unification model 
do not predict an increased obscured fraction for QSOs hosted 
by galaxies with higher Lf^. 

8. VERIFICATION OF THE DIFFERENCE BETWEEN QSO IS AND 

QS02S 

In this work, we have demonstrated that the star forma¬ 
tion properties of QSOls and QS02s are different, which im¬ 
plies a link between star-forming dust and AGN obscuration 
in quasars. To verify the observed difference between QSOls 
and QS02 is not driven by the uncertainty of photometric red- 
shifts of QS02s and the presence of starburst contaminations, 
we perform several tests on the QSOl sample and use the 
spectroscopically confirmed AGNs as a benchmark of how 
the various uncertainties might affect the observed properties 
of QS02s. 

8.1. Far-IR detection fraction 

We reiterate that the sources with a dominating starburst 
component have been removed from the final sample studied 
in this analysis based on the SED fits discussed in §3. There¬ 
fore, the higher far-IR detection fraction for QS02s is not due 
to the inclusion of starburst interlopers in the mid-IR AGN 
sample. Another possible issue that can drive the higher far- 
IR detection fraction for the QS02s is the photo-z uncertain¬ 
ties. We here examine this issue by testing the null hypothesis 
that QSOls and QS02s have exactly the same far-IR detec¬ 
tion fraction and redshift distribution. 
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Fig. 8.— The obscured fraction (i.e. the QS02 fraction) as a function of 
Ljjj. We separate far-IR detected QSOs into two bins of /.j)( with equal num¬ 
ber of sources in each bin, and plot their obscured fraction as the orange open 
squares. For the far-IR non-detected QSOs, we show the obscured fraction at 
the estimated from stacking analysis as the gray filled squares. The hori¬ 
zontal error bars show the inter-percentile range in containing 80% of the 
sample. This plot shows a monotonic increase of obscured fraction at higher 
which connects the AGN obscuration to the host galaxy star formation. 

From the AGES catalog, we select all of the sources with 
broad emission lines (the BLAGNs). In the AGES catalog 
there are a total of 1619 BLAGNs of which 181 of them 
have a SPIRE 250 pm flux greater than 20 mJy. This far- 
IR detection fraction is similar to that for the QSOl sam¬ 
ple. To test if the uncertainties in photometric redshifts can 
boost the far-IR detection fraction, we randomly scatter the 
redshifts for all BLAGNs with the conservative uncertainty 
of <r z = 0.25(1 +z) (discussed in §2.1). In each random real¬ 
ization, we select the sources with a randomly assigned red- 
shift within 0.7 < z < 1.8 and calculate their far-IR detec¬ 
tion fraction. We repeat this process 100 times and found that 
on average, the random sample would have a far-IR detec¬ 
tion fraction of only ~ 9%. This means that the uncertainty 
in the photometric redshift would actually decrease the ob¬ 
served far-IR detection fraction if the far-IR luminosity func¬ 
tion and the redshift distribution are the same for QSOls and 
QS02s. This is likely due to the negative /.'-correction _at 
250 um for galaxies with t ypical cold dust SEDs (e.g. iBlainl 
1996; Neg rello et al.1120(171) , The starburst cold dust emis¬ 
sion peaks at rest-frame ~ 100 pm which correspond to the 
SPIRE 250 pm for an object at z = 1.5. Therefore, the star- 
burst galaxies at lower redshifts are not brighter than our far- 
IR detected QSOs at observed-frame 250 pm due to the neg¬ 
ative /^-correction. High-redshift starbursts are also not likely 
to be included in our far-IR detected sample as the 250 pm 
filter probes the rapidly dropping Rayleigh-Jeans tail of the 
cold dust emission for objects beyond the redshift range of 
our sample. 

Considering the significant difference of ~ 17% between 
the far-IR detection fractions of QSOls and QS02s, we con¬ 
firm that the QS02s are indeed more likely to be hosted by 
galaxies with strong far-IR emission than QSOls. 

8.2. Far-IR luminosity 

Since the majority of QS02s have only photo-zs and lumi¬ 
nosity measurements are redshift dependent, it is of extreme 


importance to confirm that the 0.3 dex average difference 
between QSOls and QS02s is not driven by solely by the un¬ 
certainty of photo-zs . Similar to the previous subsection, we 
estimate the effect of the photo-z uncertainty by testing the 
null-hypothesis that the QS02s have the same distribu¬ 
tion as the QSOls. We again randomly assign redshifts to the 
AGES BLAGNs based on their spec-z s and the uncertainty 
upper limit o z = 0.25(1 +z). For the sources with a randomly 
assigned redshift within the 0.7 < z < 1.8 range, we perform 
the SED analysis (described in §3) using the assigned redshift 
to calculate their average . We repeat this process 100 
times and found that the average difference between the 
random samples and the original QSOl sample is 0.08 ± 0.02 
dex, which is much smaller in comparison to the observed 0.3 
dex difference between QSOls and QS02s. We also note that 
the a z = 0.25(1 +z) is a very conservative upper limit so the ef¬ 
fects of photometric redshift uncertainty are almost certainly 
smaller. The Lsf difference between QSOls and QS02s is 
not driven by the photo-z uncertainty. 


9. DISCUSSION AND SUMMARY 

Recently, a number of studies of X-ray selected AGNs have 
found that star formation does not distinguish between AGNs 
with and without obscuration for AGNs classified based on ei¬ 
ther the presence of broad emission lines dMerloni et al. |2014|) 
or the observed gas colum n density (Ah, e.g. iRovilos et all 
120121 iRosario et all 120121) . These conclusions are in direct 
contrast with our result in Fig.[6]showing a significant (> 0.25 
dex) difference in between the QSOls and the QS02s. 

To explain this apparent contradiction, we point out a fun¬ 
damental difference between these studies and the present 
work. iMerloni et all (120141) classifies the X-ray selected AGN 
into “type 1” and “type 2” sources based SED fitting or the 
detection of broad emission lines. In our R— [4.5] classifica¬ 
tion, all of the QSOls are spectroscopically selected type 1 
quasars, but it is not clear how many of our QS02s would 
have spectra with narrow emission lines similar to that of the 
optical type 2 quasars; since the large amount of dust which 
blocks most of the AGN optical continuum might very well 
block the emission lines from the narrow line region essential 
to optical AGN classifications. Also, while deep X-ray and 
optical observations have proven to be efficient at detecting 
both obscured and unobscured AGNs, it is still challenging 
to include AGNs obscured by material with column density 
reaching Nh ~ 5 x 10 24 cm -2 (i.e. Compton thick) in optical 
or X-ray selected samples. Recent works on synthetic mod¬ 
els of the cosmic X-ray backgrounds (e.g . iGilli et aD 120071 
iTreister et al.ll2009l : lBallantvne et alJl201 1 ) suggest that a sig¬ 
nificant fraction of SMBH growth might occur in a heavily 
obscured, rapidly accreting AGN phase. As we have shown 
in Fig. [7J the far-IR detected QS02s have a Lx/L^p m ratio 
consistent with AGNs that have heavily absorbed or intrinsi¬ 
cally weak X-ray emission. Therefore, the different results 
might simply be due to mid-IR selected quasars representing 
different AGN populations than the optically and X-ray se¬ 
lected AGNs. 

Indeed, mid-IR selection criteria have been shown to be 
capable of picking up sources with their nuclei obscured 
by Compton-thick material (e.g. IRovilos et all 1201 3i) . and 
moderate-luminosity mid-IR AGNs have been found to pref¬ 
erentially populate the higher-end of AGN Eddington ra¬ 
tio distributions (e.g. iHickox e t ah 2009 ), and to reside in 
more blue, star-forming host galaxies (e.g. IHickox et al.l2009i : 
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iGoulding et alJ[2014l) . In conjunction with the higher L® and 
obscuration in our mid-IR quasars, we argue that part of the 
enhancement in Ljj' of our sample when compared to X-ray 
or optical quasars is due to having different AGN populations 
based on different selection criteria and the heavy obscuration 
which possibly takes place on the scale of the host galaxy. 

In summary, we have studied mid-IR luminous quasars at 
0.7 < z < 1.8 and demonstrated a connection between obscu¬ 
ration in mid-IR selected QSOs and star formation of their 
host galaxies. This connection is supported by the higher 
AGN obscured fraction in QSO host galaxies with larger 
and the large differences in the average far-IR luminosities be¬ 
tween QSOls and QS02s. We have also shown that both the 
AGN obscuration at optical wavelengths and the AGN absorp¬ 
tion at X-ray wavelengths can be connected to the presence of 
far-IR emitting dust. These results suggest that the large scale 
gas and dust are also obscuring the central AGNs in addition 
to orientation-based small-scale obscuration, which is consis¬ 
tent with the scenario in which the rapidly growing SMBHs 


Alberts, S„ Pope, A., Brodwin, M„ et at. 2013, MNRAS, -1, 2563 
Alexander, D. M., Chary. R., Pope. A., et al. 2008. ApJ, 687, 835 
Alexander, D. M., Bauer, F. E., Brandt, W. N., et al. 2011, ApJ, 738, 44 
Antonucci, R. 1993, ARA&A, 31, 473 

Ashby, M. L. N., Stem, D., Brodwin, M., et al. 2009, ApJ, 701, 428 
Assef, R. J., Kochanek, C. S., Brodwin. M., et al. 2008, ApJ, 676, 286 
—. 2010a, ApJ, 713, 970 

Assef, R. J„ Kochanek, C. S., Ashby, M. L. N„ et al. 2010b, ApJ, 772, 26 
—. 2011, ApJ, 728, 56 

Baldwin, J. A.. Phillips, M. M„ & Terlevich, R. 1981, PASP, 93, 5 
Ballantyne, D. R., Draper, A. R., Madsen, K. K„ Rigby, J. R.. & Treister, E. 
2011, ApJ, 736, 56 

Behroozi, P. S., Wechsler, R. H., & Conroy, C. 2013, ApJ, 770, 57 
Blain, A. W. 1996, MNRAS, 283, 1340 

Brodwin, M„ Brown, M. J. I., Ashby, M. L. N., et al. 2006, ApJ, 651, 791 

Brown, M. J. I., Dey, A., Jannuzi, B. T., et al. 2007, ApJ, 654, 858 

Brusa, M., Comastri, A.. Gilli, R., et al. 2008, ApJ, 693, 36 

Burlon, D., Ajello, M., Greiner, J., et al. 2011, ApJ, 728, 58 

Canalizo, G., & Stockton, A. 2001, ApJ, 555, 719 

Ceballos, M. T„ & Barcons, X. 1996, MNRAS, 282, 493 

Chapin, E. L., Chapman, S. C., Coppin, K. E., et al. 2011, MNRAS, 411, 505 

Chary. R„ & Elbaz, D. 2001, ApJ, 556, 562 

Chen, C.-T. J., Hickox, R. C., Alberts, S., et al. 2013, ApJ, 773, 3 

Chung, S. M., Kochanek, C. S., Assef, R., et al. 2014, ApJ, 790, 54 

Dai, Y. S.. Bergeron, J., Elvis, M., et al. 2012, ApJ, 753, 33 

Dale, D. A., &"Helou, G. 2002, ApJ, 576, 159 

Del Moro, A., Alexander, D. M., Mullaney. J. R., et al. 2013, A&A, 549, 
A59 

Di Matteo, T., Springel, V., & Hemquist, L. 2005, Nature, 433, 604 
DiPompeo, M. A., Myers, A. D., Brotherton, M. S., Runnoe, J. C., & Green, 
R. F. 2014, ApJ. 787, 73 

Donley, J. L., Rieke, G. H., Rigby, J. R., & Perez-Gonzalez, P. G. 2005, ApJ, 
634, 169 

Donley, J. L., Koekemoer, A. M., Brusa, M., et al. 2012, ApJ, 748, 142 
Donoso, E., Yan, L., Stem, D., & Assef, R. J. 2014, ApJ, 789, 44 
Draine, B. T. 2003, ARA&A, 41, 46 

Drouart, G., De Breuck, C., Vernet, J., et al. 2014, A&A, 566, A53 
Eckart, M. E., McGreer, I. D., Stem, D., Harrison, F. A., & Helfand, D. J. 
2010, ApJ, 708, 584 

Eisenhardt, P. R., Stem, D., Brodwin, M., et al. 2004, ApJS, 154, 48 
Elbaz, D„ Dickinson, M„ Hwang, H. S. S., et al. 2011, A&A, 533, A119 
Fiore, F.. Puccetti, S., Brusa, M., et al. 2009, ApJ, 693, 447 
Gandhi, P, Horst, H., Smette, A., et al. 2009, A&A, 502, 457 
Georgantopoulos, I., Comastri, A., Vignali, C., et al. 2013, A&A, 555, A43 
Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Gonzalez, A. H., Brodwin, M., Brown, M. J. I., et al. 2010, American 
Astronomical Society 


are connected to dust-enshrouded starburst galaxies. 

We thank the anonymous referee for a careful reading of 
the manuscript and helpful suggestions that strengthened this 
paper greatly. We thank our colleagues on the AGES, ISS, 
SDWFS, NDWFS and the XBootes teams, and the HerMES 
team for making the data publicly available. The first Spitzer 
MIPS survey of the Bootes region was obtained using GTO 
time provided by the Spitzer Infrared Spectrograph Team (PI: 
James Houck) and by M. Rieke. We thank the collabora¬ 
tors in that work for access to the 24 p,m catalog generated 
from those data by Enteric LeFloc’h. This work was sup¬ 
ported under contract NAS8-03060. R.J.A. was supported 
by Gemini-CONICYT grant number 32120009. C.-T.J.C, 
K.N.H, and R.C.H. were partially supported by the National 
Science Foundation through grant No. 1211096. K.N.H. and 
R.C.H. were partially supported by NASA through ADAP 
award NNX12AE38G. R.C.H. acknowledges support from an 
Alfred P. Sloan Research Fellowship and Dartmouth Class of 
1962 Faculty Fellowship. C.-T.J.C was supported by a Dart¬ 
mouth Fellowship and the William H. Neukom 1964 Institute 
for Computational Science. 


Colliding, A. D„ & Alexander, D. M. 2009, MNRAS, 398, 1165 
Goulding, A. D., Alexander, D. M., Mullaney, J. R., et al. 2010, MNRAS, 
1757 “ 

—. 2011, MNRAS, 411. 1231 

Goulding, A. D.. Forman, W. R.. Hickox, R. C., et al. 2014, ApJ, 783, 40 
Griffin, M. J., Abergel, A.. Abreu, A„ et al. 2010, A&A, 518, L3 
Griffith, R. L„ & Stern, D. 2010, AJ, 140, 533 
Gmppioni, C., Pozzi, F., Rodighiero, G., et al. 2013, MNRAS, -1, 32 
Guainazzi, M., Matt, G., & Perola, G. C. 2005, A&A, 444, 119 
Hainline, K. N„ Hickox. R. C., Carroll, C. M„ et al. 2014, ApJ, 795, 124 
Harrison, C. M., Alexander, D. M., Mullaney, J. R.. et al. 2012, ApJ, 760, 
L15 

Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013. ApJ, 770, 103 
Hickox, R. C., Mullaney, J. R., Alexander. D. M., et al. 2014, ApJ, 782, 9 
Hickox. R. C., Jones, C., Forman, W. R., et al. 2007, ApJ, 671, 1365 
—. 2009, ApJ, 696, 891 

Hickox, R. C., Myers, A. D., Brodwin, M., etal. 2011, ApJ, 731, 117 
Hiner, K. D., Canalizo, G., Lacy, M., et al. 2009, ApJ, 706, 508 
Hopkins, P. F., Hemquist, L., Cox, T. J., et al. 2006, ApJS, 163, 1 
Jannuzi. B. T., & Dey, A. 1999, in ASP Conf. Ser. 191: Photometric 

Redshifts and the Detection of High Redshift Galaxies, ed. R. Weymann. 
L. Storrie-Lombardi, M. Sawicki, & R. Brunner (San Francisco: ASP), 

111 

Kennicutt, R. C. 1998, ARA&A, 36. 189 

Kenter, A., Murray, S. S., Forman, W. R., et al. 2005, ApJS, 161, 9 
Kirkpatrick, A., Pope. A., Alexander, D. M., et al. 2012, ApJ, 759, 139 
Kirkpatrick, A., Pope. A., Charmandaris, V., et al. 2013, ApJ, 763. 123 
Kochanek, C. S., Eisenstein, D. J., Cool, R. J., et al. 2012, ApJS, 200, 8 
Lacy, M., Petrie, A. O., Sajina, A., et al. 2007a, AJ, 133, 186 
Lacy, M., Ridgway, S. E., Gates, E. L., et al. 2013, ApJS, 24 
Lacy, M., Sajina, A., Petrie, A. O., et al. 2007b, ApJ, 669. L61 
Lacy, M., StorrieaARLombardi, L. J., Sajina, A., et al. 2004, ApJS. 154, 166 
Lansbury, G. B., Alexander, D. M., Moro, A. D., et al. 2014, ApJ, 785, 17 
Luo, B„ Brandt, W. N„ Alexander. D. M„ et al. 2013, ApJ, 772, 153 
Lutz, D„ Maiolino, R., Spoon, H. W. W., & Moorwood, A. F. M. 2004, 
A&A, 418, 465 

Mainieri, V., Bongiorno, A.. Merloni, A., et al. 2011. A&A, 535, A80 
Maiolino, R., Shemmer, O., Imanishi, M., et al. 2007, A&A, 468, 979 
Mateos, S., Alonso-Herrero, A., Carrera, F. J., et al. 2013, MNRAS. 434, 

941 

Mendez, A. J., Coil. A. L„ Aird, J.. et al. 2013, ApJ, 770, 40 
Merloni, A., Bongiorno, A., Brusa, M., et al. 2014, MNRAS, 437, 3550 
Mullaney, J. R.. Pannella, M„ Daddi, E„ et al. 2012, MNRAS, 419, 95 
Mullaney, J. R. R., Alexander, D. M. M., Goulding, A. D. D., & Hickox, 

R. C. C. 2011, MNRAS. 414, 1082 
Murray, S. S., Kenter, A., Forman. W. R., et al. 2005, ApJS, 161, 1 




15 


Negrello, M„ Perrotta, F.. Gonzalez, J. G.-N., et al. 2007, MNRAS, 377, 
1557 

Netzer, H. 2009, MNRAS, 399, 1907 

Netzer, H„ Lutz, D., Schweitzer, M., et al. 2007, ApJ, 666, 806 
Oliver, S. J., Bock, J„ Altieri, B„ et al. 2012, MNRAS, 424, 1614 
Page, M. J., Carrera, F. J., Stevens, J. A., Ebrero, J., & Blustin, A. J. 2011, 
MNRAS, 416, 2792 

Page, M. J., Stevens, J. A., Ivison, R. J., & Carrera, F. J. 2004, ApJ, 611, L85 
Page, M. J., Symeonidis, M., Vieira, J. D., et al. 2012, Nature, 485, 213 
Park, S. Q., Barmby, P. Willner, S. P. et al. 2010, ApJ, 717, 1181 
Polletta, M„ Tajer, M„ Maraschi, L., et al. 2007, ApJ, 663, 81 
Polletta, M. d. C., Wilkes, B. J., Siana, B„ et al. 2006, ApJ, 642, 673 
Rosario, D., Santini, P., Lutz, D., et al. 2012, A&A, 545, A45 
Rosario, D. J., Trakhtenbrot, B., Lutz, D., et al. 2013, A&A, 560, A72 
Rovilos, E., Comastri, A., Gilli, R.. et al. 2012, A&A, 546, A58 
Rovilos, E., Georgantopoulos, I., Akylas, A., et al. 2013, MNRAS, -1, 2867 
Sanders, D. B„ Soifer, B. T„ Elias, J. H„ et al. 1988, ApJ, 325, 74 
Sazonov, S., Willner, S. P., Goulding, A. D., et al. 2012, ApJ, 757, 181 
Schawinski, K., Simmons, B. D., Urry, C. M., Treister, E., & Glikman, E. 
2012, MNRAS, 425, L61 

Serjeant, S., & Hatziminaoglou, E. 2009, MNRAS, 397, 265 
Shan, H. G., & Chen, P. S. 2012, MNRAS, 421, no 


Silverman, J. D., Lamareille, F„ Maier, C., et al. 2009, ApJ, 696, 396 
Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hernquist, 
L. 2008, MNRAS, 391,481 

Stem, D„ Eisenhardt, P., Gorjian, V., et al. 2005, ApJ, 631, 163 
Stem, D., Lansbury, G. B., Assef, R. J., et al. 2014, ApJ, 794, 102 
Sturm, E., Hasinger, G., Lehmann, I., et al. 2006, ApJ, 642, 81 
Sun, A.-L., Greene, J. E., Zakamska, N. L., & Nesvadba, N. P. H. 2014, ApJ, 
790, 160 

Teng, S. H„ Brandt, W. N„ Harrison, F. A., et al. 2014, ApJ, 785, 19 
Treister, E., Urry, C. M., & Virani, S. 2009, ApJ, 696, 110 
Urry, C. M„ & Padovani, P. 1995, PASP, 107, 803 

Usman, S. M., Murray, S. S., Hickox, R. C., & Brodwin, M. 2014, ApJ, 788, 
L3 

Whitaker, K. E., van Dokkum, P. G., Brammer, G., & Franx, M. 2012, ApJ, 
754, L29 

Wilkes, B. J., Kuraszkiewicz, J., Haas, M., et al. 2013, ApJ, 773, 15 
Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 
1868 

Zakamska, N. L., Gomez, L., Strauss, M. A., & Krolik, J. H. 2008, AJ, 136, 
1607 

Zakamska, N. L„ Strauss, M. A., Krolik, J. H„ et al. 2006, AJ, 132, 1496 



